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The complex double-ridge system in the Luzon Strait in 
the South China Sea (SCS) is one of the strongest sources 
of internal tides in the oceans, associated with which are 
some of the largest amplitude internal solitary waves on 
record. An issue of debate, however, has been the specific 
nature of their generation mechanism. To provide insight, 
we present the results of a large-scale laboratory experiment 
performed at the Coriolis platform. The experiment was 
carefully designed so that the relevant dimensionless pa¬ 
rameters, which include the excursion parameter, critical¬ 
ity, Rossby and Froude numbers, closely matched the ocean 
scenario. The results advocate that a broad and coherent 
weakly-nonlinear, three-dimensional, M 2 internal tide that 
is shaped by the overall geometry of the double-ridge system 
is radiated into the South China Sea and subsequently steep¬ 
ens, as opposed to being generated by a particular feature 
or localized region within the ridge system. 


1. Introduction 

The Luzon Strait (LS), a complex double-ridge system 
between Taiwan and the Philippines, generates some of the 
strongest internal tides found anywhere in the oceans [Al¬ 
ford et aL, 2011]. Consequently, a striking characteristic of 
the South China Sea (SCS) to the west of the LS is the reg¬ 
ular observation of large amplitude Internal Solitary Waves 
(ISWs) with amplitudes up to 170 m and speeds of 3 m/s, 
taking about 33 hours to traverse the SCS basin [Alford 
et aL, 2010]. A summary of observations over a three-year 
time period reveals that these ISWs are primarily directed 
slightly north of west, with surface signatures of the ISWs 
are rarely observed east of 120.3° [Jackson, 2009]. 

In recent years, there has been a concerted effort to more 
clearly establish the early stages of the formation of ISWs in 
the SCS. The lee wave mechanism [Maxworthy, 1979], ini¬ 
tially alluded to as a candidate process [Ehhesmeyer et al, 
1991], is now considered inappropriate [Alford et al, 2010] 
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due to the lack of reported detection between the ridges, 
or observations in satellite imagery, of ISWs east of 120.3° 
[Jackson, 2009]. Instead, it is now widely accepted that 
a westward propagating, weakly-nonlinear, mode-1 internal 
tide generated at the LS steepens to form ISWs [Lien et al., 
2005; Li and Farmer, 2011]. 

A prevailing theme has been to identify generation hot 
spots within the LS. Two-dimensional models [Zhao and 
Alford, 2006; Alford et al, 2010] and numerical simula¬ 
tions employing idealized or actual transects containing pur¬ 
ported generation sites [Buijsman et al, 2010; Warn-Varnas 
et al, 2010], have been used in this regard, producing dis¬ 
cordant results. The inconsistency of the outcome of two- 
dimensional studies is most likely due to the inherently 
three-dimensional nature of the internal tide generation pro¬ 
cess, with extended sections of the two ridges playing a key 
vole [Zhang et al, 2011; Guo et al, 2011]. 

Here, we present the results of a large-scale, realistic, lab¬ 
oratory project of internal tide generation in the LS. The 
experiments, performed at the Coriolis Platform of LEGI 
(Grenoble, France), model key constituents of the genera¬ 
tion process: three-dimensional bathymetry, nonlinear strat- 
ihcation, semi-diurnal (M 2 ) tidal forcing, and background 
rotation. The goal of the experiments is to elucidate the 
geometry and characterize the nonlinearity of the radiated 
three-dimensional M 2 internal tide, thereby providing in¬ 
sight into the consequent generation of ISWs in the SCS. 
Due to the difference in the dissipation processes, laboratory 
experiments are a valuable complement to realistic three- 
dimensional numerical simulations. Any result reproduced 
by both approaches can be considered robust. We describe 
the experimental arrangement in §2 and discuss the scalings 
chosen to preserve dynamic similarity with the ocean in §3. 
Results for the M 2 barotropic and baroclinic (internal) tides 
are presented in §4 and §5 respectively. The nonlinearity of 
the internal tide is characterized in §6 before discussion and 
conclusions in §7. 

2. Experimental configuration 

A schematic of the experiment is presented in Fig. 1(a). 
It was performed at the Coriolis Platform, the world’s 
largest turntable facility for geophysical huid dynamics ex¬ 
periments, in order to incorporate the appropriate inhuence 
of background rotation and to curtail viscous effects. A 5 m 
long, 2 m wide and 0.6 m tall scale model of the LS con¬ 
structed based on topographic data provided by the Ocean 
Data Bank, National Science Council, Taiwan, was posi¬ 
tioned near the center of the tank. Partition barriers were 
introduced north and south of the model, representing Tai¬ 
wan and the Philippines, respectively. Damping material 
was afhxed to the inner walls of the tank to damp rehec- 
tions effectively, as discussed in §5. One of a pair of verti¬ 
cally oscillating prismatic plungers was used to generate the 
equivalent of a M 2 barotropic tide [Biesel and Suquet, 1951]. 
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Figure 1. (a) Schematic of the experiment. Vertical partitions (1) extend from Taiwan and The Philippines to the side 
of the tank . After hlling, the PO and SCS are separated using inserted barriers (2) and the barotropic tide is generated 
using prismatic plungers (3) located behind vertical walls (4). Dashed green lines and the dotted square indicate the 
vertical PIV planes and overhead held of view, respectively. Red stars and circles indicate the location of the CT and 
acoustic probes, respectively, (b) The experimental topography, with Taiwan (T), The Philippines (P), Itbayat Island 
(I), Sabtang & Batan Islands (S) and Babuyan Island (B) indicated, (c) Distributions of the topographic slope (black) 
and the criticality parameter (red) for the experiments (solid) and the ocean (dashed). The ocean slopes have been 
multiplied by 20 to match the laboratory range, (d) Typical winter and summer stratideations (red) in the LS, compared 
to a laboratory stratiheation (black). 


Located in the Pacihe Ocean (PO) basin behind vertical par¬ 
titions blocking internal waves emanating from the plunger 
itself, each plunger was tested individually. No differences 
were found in the resulting tidal ffows. 

A three-dimensional rendering of the experimental to¬ 
pography is presented in Eig. l(b). The horizontal and 
vertical scalings of the topography compared to the ocean 
were 1:100,000 and 1:5,000, respectively. To ensure these 
choices did not affect the dynamical processes being mod¬ 
eled, other physical quantities such as the stratiheation, tidal 
forcing and background rotation frequencies were appropri¬ 
ately rescaled, as described in §3. Only topographic features 
above 3000 m deep were reproduced, since the weak deep- 
ocean stratiheation has negligible effect on the wavelengths 
and phase speeds of the vertical modes. Several metrics were 
calculated to ascertain the quality of reproduction of the 
original ocean data set; one measure, the cross-correlation 
coefficient, was 0.90 ± 0.03, indicating a good transcription 
from the ocean to the experiment. Over the same domain, 
the distribution of slopes, V/i, presented in Eig. 1(c), con- 
hrms that the laboratory topography was indeed 20 times 
steeper than that of the ocean and with the same distribu¬ 
tion. 

A nonlinear density stratiheation was established using 
salt water introduced via a computer-controlled double¬ 
bucket system (see Eig. 1(d)). Velocity held data was ob¬ 
tained using Particle Image Velocimetry (PIV). Eor horizon¬ 
tal velocity held measurements, small particles of diameter 
350 /xm with a density in a narrow range (±1 kg/m^) cor¬ 
responding to the base of the pycnocline were used to seed 
the huid, ensuring they would clearly reveal horizontal mo¬ 
tions associated with mode-1 and mode-2 internal tides; the 
particles were illuminated from below by submerged white 
lights. Eor vertical velocity held measurements, a moveable 
laser sheet in a vertical plane cutting across the ridge il¬ 
luminated similar particles distributed uniformly in density 
to cover the entire huid depth. Local density oscillations 
in the radiated wave held were measured using conductivity 
and temperature (CT) probes, and acoustic probes on either 
side of the ridge measured the surface dehection due to the 
barotropic tide with an accuracy of 0.05 mm. 


3. Governing parameters 

There are ten dimensional parameters that characterize 
internal tide generation at the LS: the barotropic tidal fre¬ 
quency (cc) and velocity (Uo), the Coriolis frequency (/), 
the maximum and bottom value of the ocean stratiheations 
{Nm, Nb), the depth of the pycnocline ((5) and of the ocean 
(if), the vertical (ho) and across-ridge horizontal (L) scales 
of the topography, and the huid viscosity (z/), be it kinematic 
or turbulent. With two dimensions, there are thus eight gov¬ 
erning dimensionless parameters. The practical choice of 
these parameters is established [Garrett and Kunze, 2007], 
but prudence is necessary when estimating their values. In 
particular, one must differentiate between dynamics asso¬ 
ciated with global and local scale processes. We therefore 
consider both a global scale {Luzon) for the entire LS and a 
local scale (Batanes) characterized by the channel between 
Batan and Itbayat islands, which has been touted as a hot 
spot for internal tide generation. 

The eight governing dimensionless parameters are: (1,2) 
/i* = ho/H and (5* = d/H, the ratio of the topographic 
height, ho, and the pycnocline depth, (5, to the depth of the 
surrounding ocean, LT, respectively; (3) the aspect ratio of 
the topographic features ho/L; (4) the ratio of the upper 
and deep ocean stratiheations, N* = Nm/Nb] (5) the criti¬ 
cality, £* = Vh/{uj‘^ ~ P)/{^^ “ ^^)? which is the ratio 
of the local topographic slope to the internal wave ray slope; 
(6,7) the Reynolds and Rossby numbers, Re* = UoLli/ and 
Ro* = Ro//L, respectively, these being the ratios of inertia 
to viscous and Coriolis forces. Three alternative dehnitions 
of the wave nonlinearity parameter are (8): the tidal ex¬ 
cursion parameter. A* = UoluL, which is the ratio of the 
barotropic tidal excursion to the horizontal scale of the to¬ 
pography; the Long number, Lo* = Uo/Nmho, quantifying 
the tendency for upstream blocking and how separation due 
to the stratiheation; or the Eroude number, Rr* = Ro/c^, 
the ratio of the how speed to the propagation speed of 
mode-n, with n = 1 being the lowest baroclinic mode. A 
comparison of the dimensional and dimensionless parame¬ 
ters for the ocean and the lab is given in Tables 1 and 2, 
respectively. 
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Table 1. Characteristic values of key dimensional parameters for internal tide generation on the scale of the Luzon Strait {Luzon) 
and on the local scale (Batanes). The corresponding laboratory values are Labz, and Lahs^ respectively. The scaling factor relating 
laboratory values to the ocean ones is indicated at the last line. 



(jj (rad/s) 

/ (rad/s) 

Nm (rad/s) 

Afo (rad/s) 

Uq (m/s) 

5 (m) 

H (m) 

L (m) 

tiQ (m) 

n (m2/s) 

Luzon 

1.40x10-* 

5.00x10-3 

1.57x10-^ 

3.65x10-1 

0.1 

100 

3.0x103 

10^ 

1.5x103 

10-1 

Batanes 

1.40x10-* 

5.00x10-3 

1.57x10-2 

1.28x10-3 

1.0 

100 

1.5x103 

104 

8.0x102 

1 

0 

Labi, 

3.86x10-1 

1.38x10-1 

2.21 

5.15x10-2 

2.76x10-3 

0.02 

0.6 

1.0 

0.30 

10-8 

Labs 

3.86x10-1 

1.38x10-1 

2.21 

1.80x10-1 

2.76x10-2 

0.02 

0.3 

0.1 

0.16 

10-8 

Labi, 

Luzon 

2.76x103 

2.76x103 

1.41x102 

1.41x102 

2.76x10-2 

2x10-4 

2x10-4 

10-3 

5x103 

(M 

0 

1 —1 


Table 2. Characteristic values of key dimensionless parameters for the ocean and the measured experimental values. The ratio of 
laboratory to ocean values actually achieved is given at the end of the table. 



h* 

d* 

holL 

N* 


Re* 

Ro* 

A* 

Lo* 

Er* 

Er* 

Luzon 

0.50 

3.3x10-^ 

1.5x10-^ 

43.0 

[0.004 - 14] 

1.0x10^ 

2.0x10-^ 

7.14x10-^ 

2.5x10-^ 

3.6x10-^ 

7.3x10-^ 

Batan 

0.53 

6.6x10-2 

8.0x10-2 

12.3 

[0.004 - 14] 

1.0x10® 

2.0 

7.14x10-1 

2.7x10-1 

4.2x10-1 

8.3x10-1 

Labi, 

0.50 

2.7x10-2 

0.3 

42.0 

[0.01 - 14] 

2.76x103 

2.0x10-2 

7.15x10-3 

1.2x10-2 

2.1x10-2 

4.3x10-2 

Labs 

0.53 

5.5x10-2 

1.6 

1.82 

[0.01 - 14] 

2.76x103 

2.0 

7.15x10-1 

3.5x10-1 

8.9x10-1 

1.8 

Lab], 

Luzon 

1.0 

0.82 

20 

0.98 

0(1) 

2.76x10-3 

1.0 

1.00 

0.5 

0.6 

0.6 


The values in Table 2 reveals that the LS can be charac¬ 
terized as a weakly nonlinear system on the large scale, with 
a relevant influence of background rotation and topographic 
blocking; and a more nonlinear one on the local scale, with 
a large excursion parameter promoting lee waves generation 
and internal wave energy trapping by supercritical flow. The 
increase by a factor 20 in the aspect ratio and the wave ray 
slope (based on the rescaled stratification in Fig. 1(d)), pre¬ 
serving criticality as can be seen when comparing the ocean 
and laboratory distributions of e* in Fig. 1(c), is irrelevant in 
the long wave limit (equivalent to the hydrostatic approx¬ 
imation), which is expected to be valid in the generation 
process; this is of importance for the subsequent wave steep¬ 
ening leading to ISWs, however, as discussed in §6, which 
therefore is not reproduced in the experiment. Overall the 
values of the dimensionless parameters for the LS and the 



Ion 

Figure 2. M 2 tidal ellipses and horizontal velocity mag¬ 
nitude (a) in the experiment and (b) from TPXO 7.1. 
The ellipses and colorbars are scaled using the relevant 
scaling parameters listed in Table 1. (c) Vertical view 
of the velocity (arrows) and its time-averaged horizontal 
component (colormap) extracted along the W-NW tran¬ 
sect (blue dashed line) in (a), when the tidal flow from 
the SCS to the PO is maximum. 


laboratory model are consistent, except for Re^. As such, 
the nature of the oceanic internal tide generation process 
will be faithfully reproduced in the experiments. 

4. The M 2 barotropic tide 

Experiments were first performed using unstratified wa¬ 
ter to investigate solely the barotropic tide. Horizontal views 
of the flow dX z — —0.12m (corresponding to —600m in the 
ocean) are compared with TPXO 7.1 predictions [Egbert and 
Erofeeva, 2002] in Fig. 2(a) and (b); the ellipses and color- 
bars are scaled according to the parameters listed in Table 1 
so that, in principle, they should appear the same. The ori¬ 
entations of the tidal ellipses are in very good agreement; in 
both cases, the ellipses are almost flat, indicating a weak in¬ 
fluence of background rotation for most of the strait, except 
in very shallow regions where ellipses are quite large, which 
stress the importance of Coriolis effects in setting the ampli¬ 
tude of the tidal velocity near shallow features. The spatial 
distributions of the tidal amplitude are also very similar, the 
intensification of the flow with decreasing depth being well 
reproduced. Intense magnification above z = —0.12m could 
not be visualized due to the shadows cast by the laser light 
(coming from the west) striking the topography. 

Figure 2(c) presents a vertical transect of the barotropic 
tide taken along the blue dashed line in Figure 2(a). Away 
from topography a depth-uniform, nearly horizontal velocity 
is observed, as expected for barotropic tidal how. Amplih- 
cation of the how near seamounts is again visible over the 
west ridge. The acoustic surface probes revealed a phase dif¬ 
ference of ~ 155^ upon passing across the strait, which is in 
good agreement with the ~ 180^ phase shift found in TPXO 
models and numerical simulations [Niwa and Hibiya, 2004]. 
Finally, the introduction of a density stratihcation induced 
no discernible change in the depth-averaged velocities taken 
along the blue transect in Figure 2(a). Overall, we conclude 
that the barotropic tidal how responsible for internal tide 
generation in the laboratory experiments successfully repro¬ 
duced that in the ocean. 

5. The M 2 internal tide 

Figure 3(a) presents a snapshot of the east-west velocity 
in an isopycnal plane around z = —0.04m, with instanta¬ 
neous velocity vectors inside the M 2 hltered tidal ellipses. 
This view essentially contains the baroclinic velocity com¬ 
ponent by virtue of considering the time instant of how re¬ 
versal of the barotropic tide (i.e. the onset of how from the 
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SCS to the PO). Large-scale coherent waves radiate from the 
strait in a west-northwest direction. The maps of the ampli¬ 
tude and phase of the M 2 -filtered baroclinic and barotropic 
tides, which cannot be readily discriminated, are presented 
in Eig. 3(b) and (c), respectively. There is a clear modulation 
of the amplitude oriented in the ~ 285^N direction, which is 
not due to interference with internal waves reflected at the 
wall of the tank, as no signal was observed in the SCS when 
filtering solely eastward propagating waves [Mercier et a/., 
2008]. The apparent wavelength of the signal, 1.3±0.15m, 
is consistent with an internal tide dominated by the mode -1 
wavelength of 1.4 m for this stratification. A notable fea¬ 
ture in Eig. 3(c) is that a wave front, evident as a sharp line 
of constant phase, connects coherently from the central and 
southern sections of the eastern ridge to the far field west of 
both ridges, on the way smoothly adjoining the phase of the 
wave field generated by the northern and central sections 
of the western ridge; the combined wavefront coming from 
these sections of the two ridges propagates in a ~ 285^N 
direction. In the south, another distinct wavefront is ori¬ 
ented in a southwestward direction. In the northern section 
of the strait, between the two ridges, no clear phase prop¬ 
agation exists, a behavior characteristic of a standing wave 
pattern. In the absence of background rotation, although 
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Figure 3. (a) Colormap of the east-west velocity in the 
isopycnal plane at z = —0.04m at an instant of barotropic 
tide flow reversal. The local velocity direction is indicated 
by the black arrows inside the tidal ellipses, (b) Ampli¬ 
tude of the total velocity and (c) phase of the east-west 
velocity of the combined M 2 baroclinic and barotropic 
tides, filtered at the forcing frequency, (d) Same data 
as in (a) for the vertical transect indicated by the dashed 
blue line in (a), arrows indicate the in-plane velocity field. 


the wavelength of the radiated internal tide is noticeably di¬ 
minished, we found little discernible change in the direction 
of radiation, revealing the dominant role of the topography 
in shaping the radiated M 2 internal tide. 

Eigure 3(d) presents the same data as in Eig. 3(a), but 
for a vertical transect indicated by the blue dashed line. 
The structure of the velocity field takes the form of a classic 
mode -1 signal, with a vertical velocity maxima of alternat¬ 
ing signs along the length of the transect and in the vicin¬ 
ity of the pycnocline, as expected from the observations in 
Eig. 3(a). Modal decomposition confirms these results, de¬ 
termining that ~ 95% of the energy flux is in mode-1, and 
the remaining 5% in modes 2 and 3 mainly. 

6. Nonlinearity 

Although the ISWs in the SCS are best described by fully 
nonlinear models [Li and Farmer, 2011], in order to reason¬ 
ably quantify nonlinearity we assume the waves to be gov¬ 
erned by a weakly nonlinear, KdV-like equation. The non¬ 
linearity of the wave field is assessed by the ratio C/H (it is 
appropriate to scale with H since the waves are dominantly 
mode- 1 ), where C is a measure of the wave amplitude given 
by rj{x,z,t) = 77 being the vertical displacement 

of an isopycnal initially located at {x,z), and 0 the long wave 
vertical structure function of mode -1 evaluated at the initial 
depth of the isopycnal [Grimshaw et ah, 1998]. Estimations 
oiQ/H are based on time series presented in Eig. 4, extracted 
at z = —3.6 cm and coordinates ( 20 .5°N, 120 .5°E). Given 
the stratification profile, 77 is obtained from the CT probe 
recordings of the density perturbation, p in Eig. 4(a). The 
vertical velocity estimated from the vertical displacement 
{drf/dt) is compared with the vertical velocity (W) obtained 
from PIV measurements in Eig. 4(b), showing a good corre¬ 
lation between the two signals, although stronger nonlinear 
oscillations are observed in dp/dt. The east-west barotropic 
and total velocities, Fht and U respectively, are also dis¬ 
played to indicate the linear nature of the barotropic tide 
and the dominant baroclinic contribution in the total veloc¬ 
ity. Einally, the time series of in Eig. 4(c) shows that the 
wave amplitude is ~ 5% of the water height. This is consis¬ 
tent with the assumption of weak non-linearity and in good 


(a) 



15 15.5 16 16.5 17 17.5 18 18.5 19 

t/TM2 


Figure 4. Time series of (a) density perturbations p' 
from a conductivity probe located ai z = —3.6cm and co¬ 
ordinates (20.5°N,120.5°E); (b) the east-west and verti¬ 
cal velocity (U and W respectively) at the same location, 
along with the barotropic east-west velocity Vbt (exper¬ 
iments with no stratification) and the vertical velocity 
dp/dt estimated from the vertical displacements p; and 
(c) the nonlinear first internal mode amplitude ^ com¬ 
pared to the water depth H. Gray (resp. white) regions 
correspond to flow from PO to SGS (resp. SGS to PO). 
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agreement with oceanic observations and a recent analysis of 
Li and Farmer [2011] based upon the rotationally-modified 
KdV equation. 

The steepening of the radiated internal tide, producing 
sharp ISWs, is not observed in the experiments. This is 
not due to a stronger dissipation, but rather the conse¬ 
quence of the modified aspect ratio in the experiments. The 
characteristic wavelength of a soliton with amplitude Co is 
Ls ~ while the wavelength of the linear mode- 

1 is Ai ~ if/tanwith tan^ = 

Eor our experiments (Co/Lf = 0.05), the ratio Ls/Xi = 
(iL/Co)^^^ tan^ ~ 1, meaning that ISWs have the same hor¬ 
izontal scale as mode-1, while in the ocean the horizontal 
scale of ISWs is roughly 1/25 that of the mode-1 tide. The 
experiments do not preserve the scale separation between 
ISWs and the internal tide required to observe the subse¬ 
quent steepening. 

7. Conclusions and Discussion 

These experiments provide a clear demonstration that de¬ 
spite having complex three-dimensional geometry, the LS 
radiates a spatially coherent, mode-1 dominated, weakly 
nonlinear M 2 internal tide oriented in a west-northwest di¬ 
rection. The spacing of the ridges being on the order of 
the mode-1 internal wavelength for M 2 , extended sections 
of both ridges shape the westward propagating M 2 internal 
tide. We observed modulation of the amplitude of the radi¬ 
ated wave field due to interference between the barotropic 
and baroclinic tide, an effect observed in numerical models 
[Buijsman et a/., 2010]. The commensurate spacing of the 
ridges in the northern section of the LS results in a standing 
wave pattern. Being weakly nonlinear, the radiated internal 
tide is prone to steepening, and this is the origin of the ISWs 
in the SCS, although the steepening process itself is not cap¬ 
tured in these experiments. These results do not preclude 
the ability of local topographic features to produce ISWs, 
observed on a smaller scale [Ramp et a/., 2012]. We expect 
that further experimental work incorporating both the di¬ 
urnal and semi-diurnal tides will provide insight on ISWs 
with differing characteristics, such as the reported A and 
B waves [Ramp et aL, 2004]. Other potentially important 
considerations, such as the Kuroshio, could not reasonably 
be incorporated into the experiments. 
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